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Abstract 

We study the SU(3) gauge theory with Nf = 12 flavors in the fundamental 
representation by use of lattice simulations with staggered fermions. With 
a non-improved action we observe a chiral zero-temperature (bulk) transi- 
tion separating a region at weak coupling, where chiral symmetry is realized, 
from a region at strong coupling where chiral symmetry is broken. With im- 
proved actions, a more complicated pattern emerges, and in particular two 
first order transitions in the chiral limit appear. We observe that at suffi- 
ciently strong coupling the next-to-nearest neighbor terms of the improved 
lattice action are no longer irrelevant and can indeed modify the pattern ob- 
served without improvement. Baryon number conservation can be realized 
in an unusual way, allowing for an otherwise prohibited oscillating term in 
the pseudoscalar channel. We discuss the phenomenon by means of explicit 
examples borrowed from statistical mechanics. Finally, these observations 
can also be useful when simulating other strongly coupled systems on the 
lattice, such as graphene. 
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1. Introduction 



In recent years attention has been drawn to the study of conformal sym- 
metry restoration in non abelian gauge theories. On the one hand, there is 
theoretical interest in uncovering their phase diagram. On the other hand, 
the start of LHC activities creates the possibility of putting under scrutiny 
candidate scenarios for electroweak symmetry breaking, among others the 
possibility that strongly coupled dynamics govern the physics beyond the 
Standard Model. Some of these models live in a quasi-conformal region of 
the parameter space at the TeV scale, such as walking technicolor or gen- 
eralizations to composite Higgs models, or conformal symmetry might be 
thought to play a role at much higher energies. 

The main interest of these studies is of course the theory in the continuum 
limit. However, in recent years a growing amount of work has been devoted 
to the analysis of the so-called bulk transition emerging in the lattice phase 
diagram at strong bare gauge coupling, see Fig. [TJ Early studies based on 
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Figure 1: The bulk transition line in the Nf-g 2 plane of the phase diagram for SU(3) 
gauge theories with (unimproved) staggered fermions. The bulk transition separates a 
QED-like, chirally symmetric, region (S, right side) from a chirally broken phase (AS, left 
side). The data points are all for a bare lattice fermion mass of 0.025 and should ideally be 
extrapolated to the chiral limit. Data for Nf — 16 (work in progress) agree with Ref. pQ 
after mass rescaling, the point at g 2 — oo and Nf ~ 50 is from Ref. [2] and the point 
for Nf — 12 is from the present work. The end point of the bulk line is unknown. Rcfs. 
[31 HI [5] reported a further bulk transition in the chirally symmetric phase. In this paper 
we argue that next-to-nearest neighbor interactions in the improved fermion action are 
necessary for the second transition to occur. 

the strong coupling expansion of QCD predict that chiral symmetry is always 
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broken in the strong coupling limit, regardless of the number of flavors. Sup- 
port for this claim was offered by Damgaard et al. pQ, who uncovered a bulk 
transition for sixteen fundamental fermions. Later on, we have found [6] that 
QCD with eight fundamental fermions is still in the QCD phase, hence for 
Nf = 8 the chirally restoring transition has a genuine thermal - as opposed 
to bulk - nature. These findings have been further confirmed by [7j. A true 
bulk transition appears instead in QCD with twelve flavors, and the careful 
scrutiny of the region at its weak coupling side seemed consistent with exact 
chiral symmetry [HI HI- The presence of a zero temperature (bulk) transition 
between a chirally symmetric phase at weaker coupling and a chirally broken 
phase at stronger coupling is indicative of a theory being in the conformal 
window [8] . An analysis of the thermal transition in the preconformal region 
is again consistent with a critical number of flavors NJ ~ 12 |9J. A review 
of current investigations of the Nf = 12 theory, as well as a discussion of 
signatures and strategies can be found e.g. in JTUJ [TTj, [121 C3]. 

Interestingly, and amusingly, a second bulk transition was uncovered by 
us [31 Hj, between the first observed bulk transition and the weak coupling 
region, where chiral symmetry studies were carried out. It is important to 
observe that, as all our analysis was done at the weak coupling side of such 
second transition, all our conclusions on the nature of the Nf = 12 theory 
remain unaffected. Still, this was an interesting and unexpected observa- 
tion calling for further analysis. The existence of a second transition was 
confirmed by the work in j5j E] , where it was observed that the shift sym- 
metry of staggered fermions was broken in the intermediate phase. Finally, 
a recent interesting development re-examined the early strong coupling stud- 
ies: contrary to previous conclusions, it was observed that, with unimproved 
fermions, the line of bulk transitions ends for Nf ~ 51. No second transition 
was observed in this case [2]. 

A more general line of work involving quantum - or bulk - transitions 
in a particle physics environment dates back to early studies of QED at 
strong coupling. The transition in this context has been for a long time 
investigated in the hope of finding an interacting, non asymptotically free 
theory in four dimensions. Such a theory requires a second order transition 
with non-trivial exponents. Indeed, the bulk transition for QCD with a large 
number of flavors has close similarities with the QED transition and in this 
spirit, inspired by the work in [15], we have proposed [3] to search for an 
interacting UVFP at the bulk transition itself. One of us has also explored 
this possibility in the context of AdS/CFT |16j . QED-like lattice systems are 
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also being used for the simulation of strongly coupled graphene. Using an 
effective field theory description, the system can be modeled by QED in 2+1 
dimensions, whose bulk transition can be analyzed borrowing early lattice 
methods and strategies [TT] . 

In conclusion, bulk transitions are interesting for several reasons ranging 
from a diagnostic of the conformal window to fundamental QFT questions 
and the physics of condensed matter systems, such as graphene. Lattice 
methods are mandatory for studying these phenomena. And it is important 
to realize that, since we are not taking the continuum limit, lattice actions 
that are equivalent in the continuum might have substantially different fea- 
tures at finite lattice spacing. A very well known example is offered by lattice 
QED in four dimensions, where the compact and non-compact formulation 
produce a different order for the phase transition. Less investigated, but for 
many reasons interesting, is the effect of improvement. 

This note is dedicated to the study of the bulk transition for Nf = 12 
and the role of improvement. As we have anticipated, the discovery of an 
intermediate phase at strong coupling with a peculiar behavior has recently 
attracted some interest. In this note we present our results on the nature 
of this phase, and we show that it only exists when the fermion sector is 
improved. 

2. The Actions 

We simulated the SU(3) gauge theory with twelve flavors of staggered 
fermions in the fundamental representation. In order to separate the effects 
of improvement for the gauge and fermion action, we performed simulations 
for different cases, labeled A to D in Table [TJ with improvement present in the 
gauge and/or fermion sector. Many of the comparisons presented here are 
for a bare lattice fermion mass of 0.025 and a volume 16 3 x 24, as in Fig. |2} 
For some of the Actions we have explored an extended set of parameters, 
although a complete presentation of our results will appear elsewhere. 
Action A is given by: 

S = -^-Tr In M(am,U) + /3Re(l - U{V)) (1) 

where M(am, U) is the fermion matrix for the naive staggered action for a 
single flavor with mass m, (3 = Q/g 2 is the SU(3) lattice coupling and U(V) 
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Action 


Gauge Improvement 


Fermion Improvement 


A 


No 


No 


B 


Yes 


No 


C 


No 


Yes 


D 


Yes 


Yes 



Table 1: Actions used in this work: gauge improvement refers to tree-level Symanzik 
improvement in the gauge action, while fermion improvement refers to tree-level Symanzik 
improvement of the staggered fermion action, i.e. the addition of the Naik term [18 ) I19 j . 



is the trace of the ordered product of link variables along the single plaquette 
P divided by the number of colors. 

Tree- level Symanzik improvement of the gauge action leads to Action B, 

S = -^Tr In M(am, U) + £ ft(<? 2 ) £ Re(l - 17(C)) (2) 
' i=o,i eeSt 

where U(C) are the traces of the ordered product of link variables along the 
closed paths C divided by the number of colors. The Sq contains all the 
lxl plaquettes (nearest neighbors), while S\ contains all the 1x2 and 
2x1 rectangles (next-to-nearest neighbors). The couplings are defined as 
p = (5/3)/3 and ft = -(1/12)/?, where f5 = Q/g 2 is the SU(3) lattice coupling 
of the unimproved gauge action. 

Improvement of the staggered fermion action is realized according to the 
Naik prescription [TBI [19] 

S F = a 4 ^?7 M ( x )x( x )-L {ci [Upfflxix + fi) - U ] (x - fi)x(x - fi)] 

+c 2 \U^(x)U^{x + ii)U^(x + 2fi)x{x + 3/i) 
-Ul(x - n)Ul(x - 2fj)Ul(x - 3fi) X (x - 3/i)] } 

+a 4 m^2x(x)x(x) (3) 

X 

where the phase factor 7]^(x) = (— ^(^o+ai-.+s^-i) an( j action is written in 
terms of the one component staggered fermion fields x( x )- The coefficients 
Ci = 1 and C2 = reproduce the naive staggered fermion action, while the 
Naik choice C\ =9/8 and c 2 = —1/24 provides 0(a 2 ) accuracy at tree level. 
Notice that the additional Naik term involves up to third-nearest neighbor 
interactions. 
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3. Results 

The main numerical result of this work is summarized by Fig. [2j where 
the transition region for the fully improved and unimproved actions is shown. 
Two rapid crossovers are present with the improved action, while a single 
chiral symmetry breaking transition is present in the unimproved case. 
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Figure 2: The chiral condensate for the SU(3) gauge theory with Nf = 12 fundamental 
flavors as a function of 1/g 2 , with g the lattice bare coupling. We show the results for the 
unimproved action, Action A (leftmost, black) and for the improved gauge and fermion 
action, Action D (rightmost, red). Data are for am=0.025 and volume 16 3 x 24. The 
weaker coupling crossover of the improved action disappears in the unimproved case. 




3.1. Action A: the unimproved case 

Fig. [3]^a) shows the rapid crossover for the chiral condensate (left) super- 
imposed on the plaquette (right), at the bare lattice mass am = 0.025. No 
additional structure is observed in the chiral condensate. We corroborate 
these observations by showing the connected component of the chiral sus- 
ceptibility Xconn m Fig. |3](b); its behavior is as expected and no sign of an 
intermediate phase at weaker coupling and additional transitions is present. 
Given the absence of phase transitions or indications for a crossover, it is 
plausible to conclude that the weak coupling phase of this theory is continu- 
ously connected with the asymptotically free regime that admits a continuum 
limitr] If this is true, its symmetry properties are the same as the ones of 



In other words, when no phase transition occurs at the infrared fixed point (IRFP) 
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Figure 3: (a) Rapid crossover of the chiral condensate (PBP) and the plaquette for Nf = 12 
flavors with the unimproved action (Action A) as a function of the lattice coupling (3 = 6/g 2 
in the strong coupling region, for am = 0.025 and volume 16 3 x 24. (b) The connected 
susceptibility for the same parameters. 



the improved action, extensively investigated in our previous work. We then 
conclude that the rapid crossover observed for Action A in Fig. [3] should be in- 
terpreted as the finite mass remnant of a bulk chiral transition separating the 
chirally broken phase at strong coupling from the chirally symmetric phase, 
in complete analogy with the unimproved results of Ref. [1] for Nf = 16. 

3.2. Action D: the improved gauge and fermion action 

We now consider Action D, i.e. the case where both fermion and gauge 
actions are tree-level Symanzik improved. For small enough bare masses 
(am < 0.04), at the simulated volumes, two rapid crossovers are observed in 
the value of the chiral condensate: a large one at stronger coupling and a 
smaller one at weaker coupling (see Fig. [2]). Preliminary results were reported 
in [3j Hj. As expected, the transition to the chirally broken phase moves 
towards stronger couplings when the action is improved. Less expected is the 
fact that the transition appears to be realized in two steps, leading to one 
intermediate region. Fig. |4](a) shows that the crossover at stronger coupling 
becomes more pronounced as the bare mass decreases. No dependence on 
the lattice temporal extent is observed and no perturbative scaling can be 
realized [31 Hj. The mass dependence of the location of the strong coupling 
rapid crossover in Fig. [41(b) is in agreement with a linear scaling expected for 



of the theory, the strong coupling QED-like side of the IRFP should be continuously 
connected to the asymptotically free weak coupling side. 
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Figure 4: (a) Rapid crossovers in the chiral condensate (PBP) with the improved action 
as a function of the coupling /3q = 10 j g 2 in the strong coupling region for different bare 
masses, (b) The mass dependence of the critical f3 value extracted from the central point 
of the strong coupling is in agreement with a linear scaling expected for a first order 
transition. 




(a) (b) 

Figure 5: Mass dependence of the connected chiral susceptibility (left) and the discon- 
nected chiral susceptibility (right). 



a first order transition. 

The disconnected component of the chiral susceptibility shows a pronounced 
peak only in correspondence to the strong coupling rapid crossover, as shown 
in Fig. [5^b). These results indicate that the strong coupling rapid crossover 
is the one corresponding to chiral symmetry breaking. Consider now the 
crossover in the chiral condensate at weaker coupling. The hints at a jump 
become weaker as we approach the chiral limit, see Fig. [4^a) . On the other 
hand, the behavior of the chiral condensate as a function of the mass suggests 
a discontinuity in its mass derivative, which is best studied by considering 
the chiral susceptibility. 

The connected component of the chiral susceptibility exhibits near disconti- 
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Figure 6: Mass dependence of the difference of the pseudoscalar susceptibility Xtt = 
('tptp}/m and the connected scalar susceptibility. This order parameter probes J7a(1) (ef- 
fective) restoration in the continuum theory. Its behavior would suggest Ua(1) restoration 
(see caveats in main text) at the chiral transition. 



nuities at the condensate crossovers, as shown in Fig. |5[a). The magnitude of 
both discontinuities increases as the bare mass decreases. This suggests that 
the jump at weaker coupling also corresponds to a genuine phase transition 
in the chiral limit, as suggested in [3] and confirmed in [HE]. We conclude 
that we are observing two distinct phase transitions, one associated with a 
change of the slope of the chiral condensate at weaker coupling, the other 
with the chiral condensate itself at stronger coupling. 

In the continuum language the observed pattern of the susceptibilities, and 
of the continuum order parameter shown in Fig. [6j suggests Uj^(l) (effective) 
restoration at the strong coupling chiral transition. However, a proper anal- 
ysis of the axial anomaly is hampered at strong coupling by the absence of a 
conserved local flavor singlet current. It would be interesting to further ana- 
lyze the lattice non-local order parameter in the intermediate region in terms 
of the point-split staggered correlators, analogous to the finite temperature 
study in [H]. 

3.3. The spectrum 

We recall that staggered meson correlators on a lattice with temporal 
extent T and periodic boundary conditions have the general form 

C(t) = J2 A i ( e ~ mit + e- mi(T -* } ) + (-l)'i* (e-**** + e-^ (T - t} ) . (4) 

i 

For each state, the parity partner adds a component with alternating sign 
(— 1)*. This is a property of the staggered formulation and it is true for 
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all correlators with the exception of the equal mass Goldstone pseudoscalar 
correlator. For equal quark and antiquark masses, the parity partner operator 
for the Goldstone pion is proportional to a charge density operator and thus 
its vacuum expectation value is zero. 

We give an overview of our results in Fig. [7j The most salient feature in 
Fig. [7] is an oscillatory component that arises for the pseudoscalar correlator 
in the intermediate region (/3q = 10 / g 2 = 3.025). This effect was also ob- 
served by the authors of [5]. In this region chiral symmetry is exact and the 
scalar and pseudoscalar correlators should become increasingly degenerate 
by moving towards the chiral limit. 

What we see in Fig. [7J moving from weak to strong coupling (right to left), 
is as follows. In the chirally symmetric region, the pseudoscalar and scalar 
correlators are close to each other. As expected, the staggered scalar cor- 
relator has an oscillating component while the pseudoscalar has not. The 
non-horizontal shape of the ratios indicates a significant contribution from 
excited states. 

In the intermediate region, a new oscillating component arises in the 
pseudoscalar correlator, and seems to also arise in the scalar correlator for 
0o = 3.025. This is consistent with the abrupt change of slope in the mass de- 
pendence of the chiral condensate, given that the chiral susceptibility {ip?p)/m 
equals the volume integral of the pseudoscalar correlator. 

At strong coupling (/3 = 2.6) chiral symmetry is broken and the pseu- 
doscalar lightest state is the Goldstone boson of the broken symmetry, thus 
very light and largely non degenerate with the scalar state. We observe that 
the oscillating component in the pseudoscalar correlator visibly decouples. 
A second observed effect is the presence of an asymmetry under t — > T — t 
of all studied correlators in the intermediate region, i.e. Pq = 3.025. To 
highlight this asymmetry we have plotted the difference C (t) — C(T — t) for 
the pseudoscalar correlator in Fig. |8j We see that 

C(t) ^ C(T - t) for t odd 

C(t) ~ C(T - t) for t even (5) 

In other words there is a violation of staggered-time reversal symmetry. The 
asymmetry is well fitted, see Fig. [8j by the functional form 

C (1 - (-1)*) ( e - mt - e- m{T ~^) (6) 

with C ~ 1027 and m ~ 0.62 consistent with the fit of the pseudoscalar 
correlator on t «C T. 
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Figure 7: Central values of the ratios C (t) / C(t— 1) for the pseudoscalar (PS, black circles) 
and scalar (SC, red squares) two-point correlation functions for coupling values in the three 
interesting regions. From left to right: the chirally broken phase, the intermediate phase 
and the weak coupling phase. The coupling in this case corresponds (left to right) to the 
improved /3 = IO/.9 2 = 2.6, 3.025, 3.8. 
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(a) 

Figure 8: The asymmetry of the Goldstone pseudoscalar correlator in the intermediate 
region (/3 = 3.025), with superimposed the result of the fit to eq. the fitted parameters 
are m ~ 0.62 and C ~ 1027. 
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Figure 9: Two rapid crossovers are observed for the Naik improved fermion action (red), 
while one rapid crossover is observed for the Symanzik improved gauge action (black). 



One caveat is in order: it is known that such an asymmetry may typically be 
present when configurations are not thermalized or statistics is too low. For 
this reason we have increased thermalization time and statistics for this point 
to a few times the ones in the other two regions. The asymmetry persists and 
does not vary with increasing thermalization or statistics. Hence, even if the 
observed asymmetric state is a metastable state, its tunneling probability to 
the opposite asymmetry seems to be extremely low, suggesting that a seed 
is indeed stabilizing it. 

3.3.1. Disentangling the effect of fermion and gauge improvement 

In order to expose the separate effects of improvement of the fermion 
action and gauge action, we have performed two additional sets of lattice 
simulations with one improvement at a time - Action B and C in Table I. In 
Fig. [9] we show the results for the improved fermion action, or the improved 
gauge action. These results make clear that in the case under study the Naik 
improvement of the staggered fermion action is mainly responsible for the 
appearance of an intermediate region in the gauge coupling. 

4. Discussion 

The results presented above suggest that third-nearest neighbor terms in 
the Naik improved fermion action are responsible for the appearance of an 
intermediate phase. 
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This is perhaps not unexpected. It is well known that models with com- 
peting interactions may give rise to non-homogeneous structures and novel 
phase transitions. One prototypical example is the axial next-to-nearest 
neighbor Ising model, known as the ANNNI model [21]. These effects have 
not been observed at weak coupling, where non-nearest neighbor terms con- 
cur to a faster approach to the continuum limit, but might well appear at 
strong coupling when those terms become relevant. It is quite possible that 
quantitative predictions for the appearance and properties of the additional 
phase could be obtained in the framework of a strong coupling expansion 
that takes the improvement term into account explicitly - we do not pursue 
this here. 

Here, we provide a plausible argument that accounts for the appearance of 
such an intermediate phase and its peculiar properties: i) the emergence of 
an oscillating component of the staggered two-point correlation function in 
the pseudoscalar channel, and ii) the asymmetry of all two-point correlation 
functions under t — > T — t. 

The authors of [22J have considered the exactly solvable Ising chain (ID) 
of length L with next-to-nearest neighbor interactions. This example is ex- 
tremely instructive. There are two regions of parameters. In one region the 
eigenvalues of the transfer matrix are real and positive. In the other region, 
pairs of complex conjugate eigenvalues appear. Intuitively, the first region 
(region I) is where the nearest neighbor interaction is dominant, while the 
second region (region II) is where the next-to-nearest neighbor term becomes 
dominant. 

As observed in [22J, the two regions will also emerge in a Symanzik improved 
gauge action where the couplings (3 and (3i are fixed as a function of the 
inverse gauge coupling /3. In other words, it is the competition of near- 
est neighbor and next-to-nearest neighbor interactions at increasingly coarse 
lattice spacing that causes the system to enter the second region. 

The same argument can be repeated for the Naik improved staggered 
fermion action, with up to third-nearest neighbors. In this case, the emer- 
gence of complex eigenvalues of the transfer matrix can be understood by 
looking at the free lattice fermion propagator for a single flavor, given by 

SfIp)' 1 = $^7m Qsinjv - ^sin3p M j (7) 

A 1 

with — 7r/2 < p < 7r/2. The interacting theory at strong coupling can in 
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principle significantly modify the coefficient of each sine contribution. In par- 
ticular, the change of sign of the second term will induce a pair of imaginary 
poles (zero tri-momentum) in the massless dressed propagator, i.e. ghosts 
will appeal^] This signals the emergence of region II, likely the intermediate 
phase we have observed. 

It would certainly be interesting to understand more quantitatively the con- 
nection between the poles in the quark propagator - as emerging from the 
non-Hermiticity of the transfer matrix with Symanzik improvement [23] 
and the detailed structure of the two-point correlation functions in the inter- 
mediate phase. We postpone this analysis hopefully to future work. Here, 
we offer a qualitative explanation as to why a chirally symmetric phase with 
the observed exotic features can appear in a gauge theory with fermion im- 
provement. 

In general, the occurrence of an oscillatory secondary state in the pseu- 
doscalar (Goldstone) correlator with staggered fermions is forbidden by the 
baryon current conservation. With improvement of the action, the total 
fermionic current will include additional terms which in turn define a modi- 
fied form of the baryon number operator at zero chemical potential. For the 
Naik improved free fermion action this construction has been explicitly given 
by Gavai [23] . In the interacting simple construction that should suf- 

fice for our purpose starts with implementing the Kogut-Hasenfratz-Karsch 
prescription [2U 125] U(x) —> exp(/i)£7(x), W(x) — > exp(—fi)W(x) along the 
temporal direction. The total baryon number density is then 

n(ju) = d/d^logZ(n) = + n 3 (ju) (8) 

where comes from local interactions and comes from the third- 

nearest neighbor term. At vanishing chemical potential the total density 
n(fi = 0) must vanish due to baryon number conservation. This can be real- 
ized in two ways, either = 0) = n 3 (/i = 0) = 0, or rii(fi = 0) = —n 3 (fi = 
0) ^ 0. When the vanishing baryon number is realized in the second way, a 
non-zero oscillating component is allowed to appear in the (Goldstone) pseu- 
doscalar channel, as its coefficient is roughly speaking, proportional to n\. 
At the same time, n\ is also a measure of the forward-backward asymmetry. 



2 It is known that the dispersion relation for Naik improved staggered fermions always 
contains complex roots at non-zero tri-momentum. All ghosts generated by the improve- 
ment decouple in asymptotically free theories when approaching the continuum limit. 
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Hence, n\ 7^ allows an oscillating term in the pseudoscalar channel and 
a time asymmetry in all correlators; this is indeed what we observe for the 
pseudoscalar correlator and the other correlators in the intermediate region. 

Putting all the elements together, we would then arrive at this simpli- 
fied picture: with an improved staggered fermion action the occurrence of 
imaginary poles of the quark propagator (or equivalently complex eigenvalues 
of the transfer matrix) opens the possibility of intermediate phases. In the 
chirally symmetric phase, and for sufficiently weak coupling, the tendency to- 
wards degeneracy of the scalar and pseudoscalar propagators is contrasted by 
the requirement of zero baryon number; the latter forces the amplitude of the 
oscillating component in the pseudoscalar channel to vanish (0 = n ~ n\), 
while the oscillating component in the scalar channel starts appearing. To- 
wards stronger coupling, third-nearest neighbor interactions in the improved 
fermion action become increasingly relevant. Now, chiral symmetry (i.e. the 
degeneracy between scalar and pseudoscalar correlators) can still be pre- 
served by allowing n 3 = —n\ ^ 0. In this way the oscillating component 
appears in the pseudoscalar channel and the temporal asymmetry appears 
in all channels. This is the intermediate phase. When the coupling grows 
even larger, chiral symmetry is finally broken, the lightest pseudoscalar is its 
Goldstone boson, and scalar and pseudoscalar correlators can depart from 
each other. Our observations in the broken phase (see Fig. [7]) are consistent 
with a situation where the conservation of baryon number is again realized 
in the usual way. 

We note that such a characterization of the intermediate phase is consis- 
tent with the breaking of the shift symmetry discussed in [5]. In fact, the 
'partial baryon number operators' n\ and n 3 might well be directly related 
to the operators measuring the breaking of the shift symmetry. In addition, 
the presence of ghost poles in the quark propagator of the improved fermion 
action can be translated into the presence of complex eigenvalues of the im- 
proved transfer matrix, and in the case of staggered fermions the (real time) 
transfer matrix is related to the shift operator as T — Sf. 

Our study shows that the emergence of an exotic intermediate phase in the 
chirally symmetric SU(3) gauge theory with twelve fundamental (staggered) 
flavors is due to the improvement of the fermion action, which adds next- 
to-nearest neighbor interactions that compete with the local terms at strong 
coupling. 
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These observations might be of interest to model builders, when needing 
to realize exotic intermediate structures in interacting gauge models with a 
relatively simple and controlled procedure. It is also amusing to notice that 
it is possible to mimic features of a dense system and a complex action (time 
asymmetry) by working with a real action, without the sign problem. Of 
course, from the perspective of the study of the phase diagram for SU(N) 
gauge theories with many flavors, the observed features remain a peculiar 
form of lattice artifacts that should be well disentangled from the underlying 
physics of the system. 
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